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Abstract: Two very different problems that can be studied by renormalization group meth- 
ods are discussed with the aim of showing the conceptual unity that renormalization group 
has introduced in some areas of theoretical Physics. The two problems are: the ground 
state theory of a one dimensional quantum Fermi liquid and the existence of quasi peri- 
odic motions in classical mechanical systems close to integrable ones. I summarize here 
the main ideas and show that the two treatments, although completely independent of 
each other, are strikingly similar. 

§1. Introduction. 

There are few cases in which a renormalization group analysis can be performed in 
full detail and without approximations. The best known case is the hierarchical model 
theory of Wilson, [Wi70] , [WK74] . Other examples are the (Euclidean) ip 4 quantum field 
theories in two and three space-time dimensions, [W173] (for an analysis in the spirit of 
what follows see [Ga85] or [BG95]), and the universality of critical points [WF72]. In all 
such examples there is a basic difficulty to overcome: namely the samples of the fields can 
be unboundedly large: this does not destroy the method because such large values have 
extremely small probability, [Ga85] . The necessity of a different treatment of the large 
and the small field values hides, to some extent, the intrinsic simplicity and elegance of 
the approach: unnecessarily so as the end result is that one can essentially ignore (to the 
extent that it is not even mentioned in most application oriented discussions) the large 
field values and treat the renormalization problem perturbatively, as if the large fields 
were not possible. 

Here I shall discuss two (non trivial) problems in which the large field difficulties are not 
at all present, and the theory leads to a convergent perturbative solution of the problem 
(unlike the the above mentioned classical cases, in which the perturbation expansion 
cannot be analytic in the perturbation parameter). 

The problems are: 

(1) the theory of the ground state of a system of (spinless, for simplicity) fermions in 
1-dimension, [BGPS] , [BG90] , [BM95] ; 

(2) the theory of KAM tori in classical mechanics, [E196], [Ga95], [BGGM], [GGM95]. 

The two problems will be treated independently, for completeness, although it will ap- 
pear that they are closely related. Since the discussion of problem (1) is quite technical 
we summarize it at the end (in §3) in a form that shows the generality of the method 
that will then be applied to the problem (2) in §4. 

The analysis of the above examples suggests methods to study and solve several problems 
in the theory of rapidly perturbed quasi periodic unstable motion, [Ga95], [GGM99]: but 
for brevity we shall only refer to the literature for such applications. 
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§2. Fermi systems in one dimension. 



The Hamiltonian for a system of TV spinless fermions at x 1 ,...,x N enclosed in a box 
(actually an interval) of size V is: 

N ^ N 

F = E i^ A ^ A + 2A E <^ - £i) - E v ( 2A ) 

z—l i<j i—1 

where fi is the chemical potential, v is a smooth interaction pair potential, A is the 
strength of the coupling; v is a correction to the chemical potential that vanishes for 
A = and that has to be adjusted as a function of A; it is introduced in order that the 
Fermi momentum stays A-independent and equal, therefore, to pp = (2m/i)5. It is in 
fact convenient to develop the theory at fixed Fermi momentum because the latter has 
a more direct physical meaning than the chemical potential as it marks the location of 
important singukarities of the functions that describe the theory. The parameter m is 
the mass of the particles in absence of interaction. 

It is well known, [LW60] , that the ground state of the above Hamiltonian is described 
by the Schwinger functions of a fermionic theory whose fields will be denoted Vx • For 
instance the occupation number function n k which, in absence of interaction, is the simple 
characteristic function = 1 if \ k\ < pf and rik = if | k | > pp is, in general, the 
Fourier transform of S(x,t) with x = (x,t) = (x 1 — x 2 ,t\ — t?) and t = + : 

S(x) = SU^tv, x 2 ,t 2 )\ + = Jim ^ 



lti=t 2 /s^oo Tre^^ ti=tr 

(2.2) 

Formal perturbation analysis of the 2-points Schwinger function S(x) and of the n-points 
natural extensions S(x\, x 2 , ■ • • , x n ) can be done and the (heuristic) theory is very simple 
in terms of Fcynman graphs. 

The n-points Schwinger function is expressed as a power series in the couplings A, v 
^ Pl>q •5^ p ' 9 ^ {xi,x 2 , . . . , x n ) with the coefficients S^ p ' q ^ computed by considering the 
(connected) Feynman graphs composed by linking together in all possible ways the fol- 
lowing basic " 'graph elements" 

(1) p "internal 4-lines graph elements" (also called "coupling graphs") and q "internal 
2-lines graph elements" (or "chemical potential vertices" ) of the form: 




Fig. 1: The two basic building blocks ("graph elements") of the the Feynman graphs for the description 
of the ground state: the first represents the potential term (2\v) in (2.1)) and the second the chemical 
potential term (v). 

where the incoming or outgoing arrows represent ip~ or respectively, and 
(2) n single lines attached to "external" vertices x_y the first half of which oriented 
towards the vertex x and the other half of them oriented away from it: 

► — • • — < 

x x 

Fig. 2: Graphical representation of the "external" lines and vertices in Feynman graphs. 

The graphs are formed by contracting {i.e. joining) together lines with equal orientation. 
The lines emerging from different nodes are regarded as distinct: we can imagine that 
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each line carries a label distinguishing it from any other, e.g. the lines are thought to be 
numbered from 1 to 4 or from 1 to 2, depending on the structure of the graph clement 
to which they belong. So that there are many graphs giving the same contributions. 
Each graph is assigned a value which is ±(p\q\) \ p v q times a product of propagators, 
one per line. The propagator for a line joining x\ to X2 is, if x\ = (x 1 ,ti), X2 = (x 2 ,t2). 

I r e -i(k (ti-t 2 )+ k (ij-i,)) 

(2^F7 -ik + (k 2 ~p 2 F )/2m 



g(xi - x 2 ) = Xl9 ^ — 0X2 = / — - — —— 2 TT77 ^ — dkodk (2.3) 



A wavy line, see Fig. 1, joining x\ with x 2 is also given a propagator 

g{x 2 -x 1 )=v{x 2 - xJSfa-h) (2.4) 

associated with the "potential" . However the wavy lines arc necessarily internal as the 
can only arise from the first graph element in Fig. 1. 

The p + q internal node labels ( x , t) must be integrated over the volume occupied by 
the system (i.e. the whole space-time when V,j3 — > oo): the result will be called the 
"integrated value" of the graph or simply, if not ambiguous, the graph value. 

Since the value of a graph has to be integrated over the labels x = (x,t) of the internal 
nodes we shall often consider also the value of a graph $ without the propagators corre- 
sponding to the external lines but integrated with respect to the positions of all nodes 
that are not attached to external lines and we call it the "kernel? of the graph i?: the 
value of a graph d will often be denoted as Vali? and the kernel by Note that the 
kernel of a graph depends on less variables: in particular it depends only on the positions 
of the internal nodes; it also depends on the labels w of the branches external to them 
through which they are connected to the external vertices. 

Introducing the notion of kernel is useful because it makes natural to collect together 
values of graphs which contain subgraphs with the same number of lines exiting them, 
i.e. whose kernels have the same number of variables. 

The function \p^S^ is given by the sum of the values of all Feynman graphs with p 
vertices of the first type in Fig. 1 and q of the second type in Fig. 1 and, of course, n 
external vertices, integrated over the internal vertices positions. As an example consider 
the following contribution to S^ 4,2 ' : 





Fig. 3: An example of a Feynman graph: in spite of its involved structure it is far simpler than 
its numerical expression, sec (2.5). A systematic consideration of graphs as "short cuts" for formulae 
permits us to visualize more easily various quantities and makes it possible to recognize cancellations 
due to symmetries. 

The value of the above graph can be easily written in formulae: apart from a global sign 
that has to be computed by a careful examination of the order in which the -labels are 
written it is 

S^' 2 \xi,x 2 ,X8,x 9 ) = ±^y J g(xw - x 2 )g(x 3 - x 1 )g(x' 3 - x w )- 

■ g(x 4 - x' 3 ) g(x 5 - xa) g{xi - x 5 ) g(x' 5 - x' A ) g{x' A - x' s ) g(x e - x 5 )- 

(2 5) 

■g(x 4 -x 5 )g(x 5 - x 4 )g(x' 7 - x 6 ) g(x g - x' 7 ) g(x s - x 7 )S(U - t' 4 )- 

■ <K*5 - t' 5 )v(x' 3 - x 3 )S(t' 3 - t 3 )v(x' 7 - x 7 )S(t' 7 - t 7 y 

■ v(x 4 — x 4 ) v(x' b — x 5 ) dx 3 dx' 3 dx 4 dx' 4 dx$ dx' 5 dx% dx 7 dx' 7 dx\o 
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which is easily derived from the figure. And one hardly sees how this formula could be 
useful, particularly if one thinks that this is but one of a large number of possibilities 
that arise in evaluating S: not to mention what we shall get when looking at higher 
orders, i.e. at S^ p - q ^ when p is a bit larger than 2. 

Many (in fact most) of the integrals over the node variables x v will, however, diverge. 
This is a "trivial" divergence due to the fact that interaction tends to change the value 
of the chemical potential. The chemical potential is related (or can be related) to the 
Fermi field propagator singularities, and the chemical potential is changed (or may be 
changed) by the interaction: the divergences are due to the naivete of the attempt at 
expanding the functions S in a power series involving functions with singularities located 
"at the wrong places" . 

The divergences disappear if the (so far free) parameter v is chosen to depend on A as: 

oo 

v = s £ j v k \ k (2.6) 
fe=i 

with the coefficients suitably defined so that the resulting power series in the single 
parameter A has coefficients free of divergences, [LW60]. 

This leads to a power series in just one parameter A and the "only" problem left is 
therefore that of the convergence of the expansion of the Schwinger functions in powers 
of A. This is non trivial because naive estimates of the sum of all graphs contributing to 
a given order p yield bounds that grow like p\, thus giving a vanishing estimate for the 
radius of convergence. 

The idea is that there are cancellations between the values of the various graphs con- 
tributing to a given order in the power series for the Schwinger functions: and that such 
cancellations can be best exhibited by further breaking up the values of the graphs and 
by again combining them conveniently. 

The "renormalization group method" can be seen in different ways: here I am proposing 
to see it as a resummation method for (possibly divergent) power series. 

Keeping the original power series in A, u, i.e. postponing the choice of v as a function 
of A, one checks the elementary fact that the propagator g(x) can be written, setting 
k = (ko, k) e R 2 , also as: 



E E e <ap ' £ 2 h fl W(2Vs) 

h=—oo w=±l ^2 7) 

9 { i l) {k) = — r — —t—r + "ncglcgiblc corrections" 
■ - v ' -ik + uk 66 



where x (&) is a function increasing from to 1 between \pp and pp, while the functions 
X^ h '{k) are the same function scaled to have support in 2 h ~ 2 pp < \k\ < 2 h pp. This means 
that for h < it is x^(fc) = x(2 _/l fcp i ^ 1 ). The simplest choice is to take X^\k) to be 
the characteristic function of z ^pp 1 ]^ > 1 and x( z ) to be the characteristic function 
of the interval [|, 1]: 



i 


^ X 








1/2 1 





k x (1) 


1 ► 




1 



Fig. 4: A (non smooth) scaling (by a factor of 2) decomposition of unity. 
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so that 



E ^ (h) ( fc ) = 1 ( 2 - 8 ) 

h— — oo 

To avoid technical problems it would be convenient to smoothen the discontinuities in 
Fig. 4 of x and x^ 1 -* turning them into C°°-functions which in a small vicinity of the 
jump increase from to 1 or decrease from 1 to 0, this is possible while still keeping the 
scaling decomposition (2.8) (i.e. with X^(fc) = x(k))- However the formalism that this 
smoothing would require is rather havy and hides the stucture of the approach; therefore 
we shall continue with the decomposition of unity in (2.8) with the sharply discontinuous 
functions in Fig. 4, warning (c.f.r. footnote 2 below) the reader when this should cause 
a problem. 

The "negligible terms" in (2.7) are terms of a similar form but which are smaller by 
a factor 2 h at least: their presence does not alter the analysis other than notationally. 
We shall henceforth set them equal to because taking them into account only introduces 
notational complications. 

The above is an infrared scale decomposition of the propagator g(x): in fact the prop- 
agator g( h ) contains only momenta k of 0(2 h pp) for ft < while the propagator g^> 
contains all (and only) large momenta (i.e. the ultraviolet part of the propagator g(x)). 
The representation (2.7) is called a quasi particles representation of the propagator and 
the quantities lo pf are called a quasi particles momenta. The function gffi is the "quasi- 
particle propagator on scale h" . 

After extracting the exponentials e l — PF - from the propagators the Fourier transforms 
<jul\k) of g^(x) will no longer be oscillating on the scale of pf and the variable k will 
have the interpretation of "momentum measured from the Fermi surface" . 

The mentioned divergences are still present because we do not yet relate A and v. they 
will be eliminated temporarily by introducing an infrared cut-off. i. e. by truncating the 
sum in (2.7) to h > —R. We then proceed keeping in mind that we must get results 
which are uniform as R — > oo: this will be eventually possible only if v is suitably fixed 
as a function of A. 

Writing g(x) = Z^g^(x) + Z^g^(x) with Z x 1 and g^ m 1 being defined in 
general, see (2.7), as: 

m 

g {<™)( x )= E ^ h e l - PF - g [ ^(2 h p F x), m<0 (2.9) 

h=-R w=±l 

each graph can now be decomposed as a sum of graphs each of which with internal lines 
carrying extra labels "1" and "w" or "< 0" and "w" (signifying that the value of the 
graph has to be computed by using the propagator Z 1 ~ 1 g^\x — x') or Zi^ 1 g^°\x — x') 
for the line in question, if it goes from x' to x). 

We now define clusters of scale 1: a "cluster" on scale 1 will be any set C of vertices 
connected by lines bearing the scale label f and which are maximal in size (i.e. they are 
not part of larger clusters of the same type) . Wavy lines are regarded as bearing a scale 
label 1. The graph is thus decomposed into smaller graphs formed by the clusters and 
connected by lines of scale < 0: it is convenient to visualize the clusters as enclosed into 
contours that include the vertices of each cluster as well as all the lines that connect 
two vertices of the same cluster. The latter can be naturally called lines internal to the 
cluster C. 

The integrated value of a graph will be represented, up to a sign which can be determined 
as described above, as a sum over the quasi particles labels to of the cluster lines and 
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as an integral over the locations of the inner vertices of the various clusters lines. The 
integrand is a product between 

(a) the kernels associated with the clusters d and depending only on the locations 
of the vertices inside the cluster Ci which are extremes of lines external to the cluster 
and on the quasi particles labels u of the lines that emerge from it, 1 and 

(b) the propagators 1 g^^ corresponding to the lines that are external to the clusters 
(in the sense that they have at least one vertex not indide the cluster). 

We now look at the clusters C that have just \C\ = 2 or |C| = 4 external lines and that 
are therefore associated with kernels Kc{{xj, w j}j=i,2; C) or Kc({xi, <±l i \i=i,...A\C). 
Such kernels, by the structure of the propagators, see (2.7) and (2.9), will have the form: 

K C = e ipFi ^ ^ Xi) K c {{ Xj , u j b=i >2 ,...,|c|) (2-10) 

where Xj are vertices of the cluster C to which the entering and exiting lines are attached; 
the cluster may contain more vertices than just the ones to which the external lines 
are attached: the positions of such "extra" vertices must be considered as integration 
variables (and as integrated) , and a sum is understood to act over all the quasi particles 
labels of the internal lines (consistent with the values of the external lines ojj's). 
If \C\ = 2,4 we write the Fourier transform at k = (fco, k) of the kernels Kq{- ■ ■)'■ 

Z 1 2- 1 v£ ) 5^ i ^ 2 + Z 1 (-ih $ ) + ^fc4 1) )<W- 2 + "remainder" 
z i^ { C 5 ^ 1 + ^ 2 + ^ 3 + ^ i =o + "remainder" 

where the first equation (\C\ — 2) is a function of one k only while the second equation 
(|C = 4) depends on four momenta k: one says that the remainders are obtained by 
"subtracting from the kernels their values at the Fermi surface" or by collecting terms 
tha do not conserve the quasi particles momenta (like terms with S u , 2 .- U j 2 in the first 
equation or with lo 1 + ...lj 4 ^0ui the second). 

The remainder contains various terms which do not have the form of the terms explicitly 
written in (2.11): a form which could be as simple as uj ■ kS^ ^ but that will in 
general be far more involved. 

In evaluating graphs we imagine, as described, them as made with clusters and that the 
graph value is obtained by integrating the product of the values of the kernels associated 
with the graph times the product of the propagators of the lines that connect different 
clusters. 

Furthermore we imagine to attach to each cluster with 2 external lines a label indicating 
that it contributes to the graph value with the first term in the decomposition in (2.11) 
only (which is the term proportional to v^ 1 " 1 ), or with the second term (which is propor- 
tional to (— ik^C,^ + uj_ka^) or with the remainder. This is easily taken into account 
by attaching to the cluster an extra label 1,2 or r. 

Likewise we imagine to attach to each cluster with 4 external lines a label indicating 
that it contributes to the value with the first term in the decomposition in (2.11) only 
(which is the term proportional to A*- 1 -*), or with the remainder. This is again easily 
taken into account by attaching to the cluster an extra label 1, r. 

The label r stands for "remainder term" or "irrelevant term" , however irrelevant does not 
mean neglegible, as usual in the renormalization group nomenclature, (on the contrary 
they are in a way the most important terms) . 

The next idea is to collect together all graphs with the same clusters structure, i.e. which 
become identical once the clusters with 2 or 4 external lines are "shrunk" to points. Since 



By definition the kernel Kq also involves integration over the locations of its inner vertices and the 
sum over the quasi particle momenta of the inner propagators. 
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the internal structure of such graphs is different this means that we are collecting together 
graphs of different perturbative order. 

In this way we obtain a representation of the Schwinger functions that is no longer a 
power series representation and the evaluation rules for graphs in which single vertex 
subgraphs (or single node subgraphs) with 2 or 4 external lines have a new meaning. 
Namely a four external lines vertex will mean a quantity ZfX' equal to the sum of Z\\% ] 
of all the values of the clusters C with 4 external lines and with label 1. 

The 2 external lines nodes will mean 

e i{u^-^x^ F ^^Z x (v' + C'd t - ia'wd s )8{x l - x 2 )) (2.12) 

where again v 1 or (',a' are the sum of the contributions from all the graphs with 2 
external lines and with label or 1 respectively. 

It is convenient to define 6' = (' — a' and to rewrite the 2 external lines nodes contri- 
butions to the product generating the value of a graph simply as: 

^(^irM.i*^^ ^ 2 Z 1 (2iy' + 5'd t + a'(d t - u ■d^))5(x 1 -x 2 ) (2.13) 

One then notes that this can be represented graphically by saying that 2 external lines 
nodes in graphs which do not carry the label r can contribute in 3 different ways to the 
product determining the graph value. The 3 ways can be distinguished by a label 0, 1' 
and z corresponding to the three addends in (2.13). 

Any graph without z-type of nodes can be turned into a graph which contains an 
arbitrary number of them, on each line connecting the clusters. And this amounts to 
saying that we can compute the series by imposing that there is not even a single vertex 
with two external lines and with label z simply by modifying the propagators of the lines 
connecting the graphs: changing them from Z-j -1 ^- ) to Z^ 1 ^- ) with 

Z = Z 1 (l + a') (2.14) 

This can be seen either elementarily by remarking that adding values of graphs which con- 
tains chains of nodes with label z amounts to summing a geometric series (i.e. precisely 
the series X)fcLo( — l) fc (a') fe = (1 + a') -1 or , much more easily, by recalling that the 
graphs are generated by a formal functional integral over Grassmanian variables and 
checking (2.14) from this remark without any real calculation, see [BG95]. In the first 
approach care is needed to get the correct relation (2.14) and it is wise to check it first in 
a few simple cases (starting with the "linear" graphs which only contain nodes with one 
entering line and one exiting line, see Fig. 1: the risk is to get Z = Z\ (1 — a 1 ) instead 
of (2.14)). 2 

Correspondingly we set: 



It is at this point that using the sharply discontinuous ^-functions would cause a problem. In fact if 
one uses the smooth decomposition (2.14) is no longer correct: namely it would become Zq = Z\ (1 + 
X (k)a) with the consequence that Z\ would no longer be a constant. At this point there are 
two possible ways out: the first is to live with a Zq which dpends on k and with the fact that the 
quantities introduced below Zj, j < —1, will also be fe-dependent; this is possible but it is perhaps 
too different from what one is used to in the phcnomcnological rcnormalization group approaches in 
which quantities like Zj arc usually constants. The other possibility is to modify the propagator on 

scale from Z~V- 0) (£) to ( k ) (1 + a')/(l + x (0) ( k ) a'). The second choice implies that 

will no longer be exactly x^ih.) / (—ikQ + w ■ k) but it will be gradually modified as h decreases 
and the modification has to be computed step by step. This is also unusual in the phcnomcnological 
rcnormalization group approaches: the reason being simply that in such approaches the decomposition 
with sharp discontinuities is always used. The latter is not really convenient if one wants to make 
estimates of large order graphs. Here this will not be a poblem for us because we shall not do the 
technical work of deriving estimates. In [BG90] as well as in [BGPS] the second choice has been 
adopted. 
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2.15 



V 




A<°> = |a' 



(2.15) 



We can now iterate the analysis: we imagine writing the propagators of the lines con- 
necting the clusters so far considered and that we shall call clusters of scale 1 as: 

= ^ (0) + (2-16) 

and proceed to decompose all the propagators of lines outside the clusters of scale 1 into 
propagators of scale or of scale < — 1. 

In this way, imagining all clusters of scale as points, we build a new level of clusters 
(whose vertices are either vertices or clusters of scale 0): they consist of maximal sets of 
clusters of scale 1 connected via paths of lines of scale 0. 

Proceeding in the same way as in the above "step 1" we represent the Schwinger func- 
tions as sums of graph values of graphs built with clusters of scale connected by lines 
with propagators on scale < —1 given by Zq 1 g^- -1 ^ and with the clusters carrying labels 
1,2 or r. Again we rearrange the 2-external lines clusters with labels 1,2 as in (2.13) 
introducing the parameters 5', A', a', v' and graphs with nodes of type z by defining a 1 
in an analogous way as the previous quantity with the same name (relative to the scale 
1 analysis). 

The one vertex nodes of such graphs with 2 or 4 external lines of scale < —1 will 
contribute to the product defining the graph value, a factor Z-\lv^~ x ^ or Z-\5 << ~ v> or 
Z 2 _ 1 \^ 1 ^ > (while the propagators in the clusters of scale 1 and the 2 or 4 nodes with two 
lines of scale emerging from them retain the previous meaning). Again one sets: 



and correspondingly we set: 

^ = 2^V, = ^S>, At" 1 ) = |±A' (2.18) 

Zjq Zjq 

and now we shall only have graphs with 2 or 4 external lines clusters which carry a label 
0, 1' or r as in the previous analysis of the scale 1 and the propagators connecting clusters 
of scale changed from Z^ 1 ^- -1 ) to ZZ\g t "~^ ■ 

Having completed the step we then "proceed in the same way" and perform "step —1" 
and so on. 

One can wander why the choice of the scaling factor 2 in (2.13) and (2.18) multiplying 
the ratio of the renormalization factors in the definition of the new Vj or, for that matter, 
why the choice of 1 for the definition of the new 5j , Xj : these are dimensional factors 
that come out naturally and any attempt at modifying the above choices leads to a beta 
function, defined below, which is not uniformly bounded as we remove the infrared cut- 
off. In other words: different scalings can be considered but there is only one which is 
useful. It could also be found by using arbitrary scaling and then look for which one the 
estimates needed to get a convergent expansion can be made. 

The conclusion is a complete rearrangement of the perturbation expansion which is now 
expressed in terms of graphs which bear various labels and, most important, contain 
propagators that bear a scale index which gives us information on the scale on which 
they are sizably different from 0. The procedure, apart from convergence problems, leads 
us to define recursively a sequence A^-* , 6^ , , Zj of constants each of which is a sum of 
a formal power series involving values of graphs with 2 or 4 external lines. The quantities 
g — (A^ , v^) can be called the running coupling constants while Zj can be called 
the running wave function renormalization constants: here j = 1,0, —1, —2, . . .. 



2.17 



Z_i = Z {l + a') 



(2.17) 
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Of course all the above is nothing but algebra, made simple by the graphical represen- 
tation of the objects that we wish to compute. The reason why it is of any interest is 
that, since the construction is recursive, one derives expressions of the g ,Zj in terms 
of the g , Z n with n > j: 



Z 



1 + B 'j(S. j+v 9 j ,---,9 ) 



9 j = ^9_ j+1 +C j {g_. +v g_.,...,g_ Q ) 



(2.19) 



where Aj is a matrix: 




(2.20) 

and the functions B'^ C_j are given by power series, so far formal, in the running cou- 
plings. The expression of Zj+i/Zj can be used to eliminate such ratios in the second 
relation of (2.19) which therefore becomes 

Zt±± = l + B' j (g g g ) 

9 . = A o . M + B Ag . , , , o.,..., o.) 

2-j — J + l — ■> y — 3 + 1 —3 '2.0' 

where A is the diagonal matrix with diagonal (1,1,2). The scalar functions B'j and the 
three components vector functions i? ■ = (Bj_i, Bj^, Bj >3 ) are called the beta functional 
of the problem. 

There are two key points, which are nontrivial at least if compared to the above simple 
algebra and which we state as propositions 

Proposition 1 (regularity and boundedness of the beta function): Suppose that there is 
e > such that g . | < e, \Zj/Zj-\ — 1| < e for all j < 1 then if e is small enough the 

power series defining the beta functionals converge. Furthermore the functions B 'j , B'^ are 
uniformly bounded and have a dependence on the arguments with label j + n exponentially 
decaying as n grows, namely there exist constants D, k such that if G = (<7 • • • , 9 ) 

and G' = (g'- +1 , ■ ■ ■ , g' Q ) with G and F_ differing only by the (j + n)-th "component" 
d = q' — a . , 7^0, then for all j < and all n > 

— S-j+n —3+n J — — 

\Bj{G)\, |S;(G)| <De 2 

\B 3 (G')-B 3 (G)\, |^(G')-S;(G)|< J De— e\d\ 

if e is small enough: i.e. the "memory" of the "beta functionals" Bj,Bj is short ranged. 
The Schwinger functions are expressed as convergent power series in g in the same 

domain \ g \ < s. 

The difficult part of the proof of the above proposition is to get the convergence of the 
series under the hypotheses \ g \ < s, \Zj/Zj_i — 1| < e for all j: this is possible because 
the system is a fermionic system and one can collect the contributions of all graphs of a 
given order k into a few, i.e. not more than an exponential in k, groups each of which 
gives a contribution that is expressed as a determinant which can be estimated without 
really expanding it into products of matrix elements (which would lead to bounding the 
order A; by a quantity growing with k\) by making use of the Gram-Hadamard inequality. 
Thus the that is in the definition of the values compensates the number of labels 
that one can put on the trees and the number of Feynman graphs that is also of order 
k\ is controlled by their representability as determinants that can be bounded without 
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generating a k\ via the Hadamard inequality. The basic technique for achieving these 
bounds is well established after the work [Le87]. A second non trivial result is 

Proposition 2 (short range and asymptotics of the beta function): 

Let G° = ( g , g , . . . , g) with g = (A, 8, v) then the function B_j(G°) defines an analytic 

function of g , that we shall call "beta functional", by setting 

(3(g) = lim B_j(G a ) (2.23) 

— — J^-oo 

for \g \ < e. The limit is reached exponentially \ (3 (g) — B_-(G )\ < e 2 De~ K ^, for some 
k > 0, D > provided \g \ < e. 

Finally the key result, [BG90], [BGPS], is 

Proposition 3 (vanishing of the beta function): If g — (A, 8, 0) then the functions 
P ( 9 ) = provided \g \ < e. Furthermore for some D, k > it is, for all j < 

B *(9 j+V ■ ■ ■ > S ) = ^+1^+1^3(5 ■ ■ ■ , £ ) + e"*B? 3 (g j+v . . . , g Q ) 
l^s(£ J+1 .-,£ )l<A ±,< ■ » 2 

provided, for h = 0, . . . , j + 1, \ g | < e. 

The above propositions are proved in [BG90], [BGPS], [BM00]. The vanishing of 
(3 (A, 8, 0) is proved in a rather indirect way. We proved that the function (3 is the 
same for the model (2.1) and for a similar model, the Luttinger model, which is exactly 
soluble; but which can be also studied with the technique described above: and the only 
way the exactly soluble model results could hold is to have (3 = . 

The vanishing of the beta function seems to be a kind of Ward identity: it is easy to 
prove it directly if one is willing to accept a formal proof. This was pointed out, after 
the work [BG90], in other papers and it was believed to be true probably much earlier 
in some equivalent form, see [So79] ; note that the notion of the beta function is intrinsic 
to the formalism of the renormalization group and therefore a precise conjecture on it 
could not even be stated before the '70s; but of course the existence and importance of 
infinitely many identities had already been noted. 

Given the above propositions one shows that "things go as if" the recursion relation for 
the running couplings was, up to exponentially small corrections, a simple memoryless 
evolution g ■_ 1 = P(g ■) + 0(e~ K ^)\ the propositions say in a precise way that this 
is asymptotically, as j — > — oo, true. This tells us that the running couplings \j, Sj 
stay constant (because (3\,(3<i vanish): however they in fact tend to a limit as j —> — oo 
exponentially fast because of the corrections in the above propositions, provided we can 
guarantee that also i/j J -^_ oo > exponentially fast and that the limits of \,,5j do not 
exceed e (so that the beta functionals and the beta function still make sense). 

It is now important to recall that we can adjust the initial value of the chemical 
potential. 3 This freedom corresponds to the possibility of changing the chemical po- 
tential 'correction" v in (2.1) and tuning its value so that Vh — ► as h — ► — oo. 

Informally if v a is chosen "too positive" then Vj will grow (exponentially) in the positive 
direction (becoming larger than e, a value beyond which the series that we are using 
become meaningless); if v is chosen "too negative" the Vj also will grow (exponentially) 



Which is a "relevant operator" , in the sense that if regarded as a running coupling it is roughly multiplied 
by 2 at each change of scale, i.e. "j—i ~ 
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in the negative direction: so there is a unique choice such that Vj can stay small (and, 
actually, it can be shown to converge to 0. 4 

The vanishing of the beta function gives us the existence of a sequence of running 
couplings g — (Xj,8j,i/j) which converge exponentially fast to (A_oo, (5_oo, 0) as j — > 
-co if v are conveniently chosen: and one can prove that X_ 00 ,8- 00 ,i' a are analytic in 
A for A small enough, [BGPS]. 

In this way one gets a convergent expansion of the Schwingcr functions: which leads to 
an essentially complete theory of the one dimensional Fermi gas with spin zero and short 
range interaction. 

§3. The conceptual scheme of the renormalization group approach followed 
above. 

The above schematic exposition of the method is a typical example of how one tries 
to apply the multiscale analysis that is commonly called a "renormalization group ap- 
proach" : 

(1) one has series that are easily shown to be finite order by order possibly provided 
that some free parameters are suitably chosen ( "formal renormalizability theory" : this 
is the proof in [LW60] that if Vh in (2.2) are suitably chosen we obtain a well defined 
perturbation series in powers of A. 

(2) However the series even when finite term by term come with poor bounds which 
grow at order k as k\ which, nevertheless are often non trivial to obtain (although this 
not so in the case (2.1) discussed here unlike the case discussed in the next sections). 

(3) One then tries to reorganize the series by leaving the original parameters (A, v) in 
the present case as \i is fixed) as independent parameters and collecting terms together. 
The aim being to show that they become very convergent power series in a sequence of 
new parameters, the "running couplings" ^ h \S^ and A^*) in the present case, under 
the assumption that such parameters are small (they are functions, possibly singular, of 
the initial parameters of the theory, A, v in the case (2.1), as 8 has to be imagined to be 
0). 

(4) The running couplings, essentially by construction, also verify a recursion relation 
that makes sense again under the assumption that the parameters are small. This relation 
allows us to express (if it makes sense) successively the running couplings in terms of the 
preceding ones: the running couplings are ordered into a sequence by "scale labels" h = 
1,0, —2, . . .. The recursion relation is interpreted as an evolution equation for a dynamical 
system (a map defined by the beta function(al)): it generates a "renormalization group 
trajectory" (the sequence (Xh,8h, Vh) out of the original parameters A, v present in (2.1), 
as 8 has to be taken as 0). 

(5) One then shows that if the free parameters in the problem, (i.e. X,v in (2.1)) are 
conveniently chosen, then the recursion relation implies that the trajectory stays bounded 
and small, thus giving a precise meaning to (2.2)), and actually the limit relation holds 
(\( h \ S^ h \ v^) h ^_ 00 ' (^oo, Soo, 0) (this is achieved in the above Fermionic problem by 
fixing v as a suitable function of A, sec [BGPS]. 

(6) Hence the whole scheme is self-consistent and it remains to check that the expressions 
that one thus attributes to the sum of the series are indeed solutions of the problem that 
has generated them: not unexpectedly this is the easy part of the work, because we have 
always worked with formal solutions which "only missed, perhaps, to be convergent" . 

(7) The first step, i.e. going to scale is different from the others as the propagators 
have no ultraviolet cut off (see the graph of in Fig. 2). Although there are no 



A simplified analysis is obtained by "neglecting memory corrections" i.e. using as a recursion relation 
g. = ^9 j+1 + £(£j +1 ) with 0(9_) verifying (2.24): this gives that Xj , 8j J _,_ OQ ' (^-00 . s -oo) ex- 
ponentially fast and va —■ — 3 — >0 exponentially fast provided uq is suitably chosen in terms of Xq.Sq: 
otherwise everything diverges. 
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ultraviolet divergences the control of this first step offers surprising difficulties (due to 
the fact that in the direction of ko the decay of the propagators is slow making various 
integrals improperly convergent): the analysis is done in [BGPS] and [GS93]. 
Note that the above scheme leaves room for the possibility that the running couplings 
rather than being analytic functions of a few of the initial free parameters are singular: 
this does not happen in the above fermionic problem because some components of the 
beta function vanish identically: this is however a peculiarity of the fermionic models. In 
other applications to field theory, and particularly in the very first example of the method 
which is the hierarchical model of Wilson, this is by far not the case and the perturbation 
series are not analytic in the running couplings but just asymptotic in the actual free 
parameters of the theory. The method however "reduces" the perturbation analysis to a 
recursion relation in small dimension (namely 3 in the case (2.1)) which is also usually 
easy to treat hcuristically. The d = 2 ground state fermionic problem (i.e. (2.1) in 2 
space dimensions) provides, however, an example in which even the heuristic analysis is 
not easy. 

In the following section we discuss another problem where the beta function does not 
vanish, but one can guarantee the existence of a bounded and small solution for the 
running couplings thanks to a "gauge symmetry" of the problem. This is an interesting 
case as the theory has no free parameters so that it would not be possible to play on them 
to find a bounded trajectory for the renormalization group running constants. This also 
illustrates another very important mechanism that can save the method in case there 
seemed to be no hope for its use, namely a symmetry that magically eliminates all terms 
that one would fear to produce "divergences" in formal expansions. Again the case 
studied is far from the complexity of gauge field theory because it again leads to the 
result that the perturbation series itself is summable (unlike gauge field theories which 
can only yield asymptotic convergence) : but it has the advantage of being a recognized 
difficult problem and therefore is a nice illustration of the role of symmetries in the 
resummation of (possibly) divergent series and the power of the renormalization group 
approach in dealing with complex problems. 

§4. The KAM problem. 

Consider d rotators with angular momentum A = (A\, . . . , Ad) G R d and positions 
a = (ai, . . . , ad) G T d = [0, 2ir] d ; let J > be their inertia moment and suppose that 
ef(a) is the potential energy in the configuration a, which we suppose to be an even 
trigonometric polynomial (for simplicity) of degree N. Then the system is Hamiltonian 
with Hamiltonian function 

n = ^jA 2 + sf(a) (4.1) 
giving rise to a model called "Thirring model". 5 



(1) The global canonical transformations C of ij" x with generating functions S(A, a) = N A 1 ■ a + 
7 (a )■ A'+ip( a) parameterized by an integer components non singular matrix JV, and analytic functions 
ff(a),/(a) leave invariant the class of Hamiltonians of the form H = (A.M(a)A) /2+ A ■ g (a )+/(«)• 
The subgroup CL^(R) of the global canonical coordinate transformations C was (remarked and) used 
by Thirring so that (4.1) is called the "Thirring model", see [Th83]. 

(2) The function H_ e {ip) in (4.2) must have zero average over j/> or, if V — * i> +5£o*> over time: hence 
the surviving quasi periodic motions can be parameterized by their spectrum uj_ or, equivalently, by 
their average action A_ = Ju .. The "spectral dispersion relation" between the average action A Q 
and the frequency spectrum is not twisted by the perturbation. Furthermore the function E(t£ ,J) can 

be taken monotonically increasing J: J~ 1 is called the "twist rate. The latter two properties motivated 
the name of "twistless motions" given to the quasi periodic motions of the form (4.2) for Hamiltonians 
like (4.1). 

(3) The invariance under the group CL d (R) has been used widely in the numerical studies of the best 
trcshold value e( uj , J) and a deeper analysis of this group would be desirable, particularly a theory of 
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For e — motions are quasi periodic (being t — > (A , a + w f) with cj = J _1 ;Aq) 
and their "spectrum" lj fills the set So = -R d of a/Z vectors w : there is a 1-to-l corre- 
spondence between the spectra lo and the angular momenta A . 

Question: If e ^ can we /ind, given lu (zSq a perturbed motion, i.e. a solution of 
the Hamilton equations of (4-1), which has spectrum oj and, as £ — > 0, reduces with 
continuity to the unperturbed motion with the same spectrum? or less formally: which 
among the possible spectra lo G So survives perturbation? 

Analytically this means asking whether two functions H Ail)), h £ (ip) on T d exist, are 
divisible by £ and are such that if A = Ju and if we set 

A =A Q + H Aip) 

~ -° ~ eK ^\ yj ET d (4.2) 
a =rp_ + h e (ij)) ~ 

then ip — ► ip + LJ t yields a solution of the equations of motion for e small enough. 
It is well known that in general only "non resonant" spectra can survive: for instance 
(KAM theorem) those which verify, for some r, C > 

1^0 ' < C| ^ r <E Z d = {integer vectors}, v^O (4.3) 

for some C, t > (Diophantine vectors) and we restrict for simplicity to such vectors. 
Furthermore given uj_ the perturbation size £ has to be small enough: |e| < £o(lj , J). 
To find H_ e , h e we should solve the equation (setting J = 1), 

ioLo ■ d ± f h(±) = -e{9*f)i± + hit)) ( 4 - 4 ) 

and if such h is given then, setting /i £ (V0 = h(ip), H-eii*) = (^o ' 9) h(ip) onc 
checks that (4.2) has the wanted property (i.e. ip — > ip + co t is a motion for (4.1)). 

To an exercized eye (4.4) defines the expectation value of the 1-particle Schwinger 
function of the euclidean field theory on the torus T d for two vector fields F ± ( ip ) = 
(F^ 1 (tp ),..., (ip)) whose free propagator is 

H1> -■>!>') 



and the interaction Lagrangian is, see [Ga95] 



{F nt)F^it'))=^^ m -— £ . J (4-5) 

^ ' „-/-c\ \ — o — ' 



C{F)=e( F + (iP)- d a (iP + F~(iP))diP 

JT d ~ ~ ~ ~ ~ 



(4.6) 



If Po(dF_) the "functional integral" defined by Wick's rule with propagator (4.5) then 

fF-(iP)e s f T ^ + ^ d -^ + ^»P (dF) 
h(ip) = =7 4.7 

J e eJ Td ^( ± y0^f( t+ F- (±))po{dE) 

At first sight this is a "sick field theory" . Not only the fields F ± ( ip ) do not come from 
a positive definite propagator, hence (4.7) has to be understood as generating a formal 
expansion in e of h with integrals over F_ being defined by the Wick rule, but also the 
theory is non polynomial and naively non renormalizable. 

The "only" simplification is that the Feynman diagrams of (4.7) are (clearly) tree- 
graphs, i.e. looplcss: this greatly simplifies the theory which, however, remains non 
renormalizable and non trivial (being equivalent to a non trivial problem). 



its unitary representations. 
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It is not difficult to work out the Feynman rules for the diagrams expressing the fc-th 
order coefficient of the power series expansion in e of (4.7). 

Consider a rooted tree with k nodes: the branches are considered oriented towards the 
root which is supposed to be reached by a single branch and which is not regarded as a 
node of the tree (hence the number of nodes and the number of branches are equal). 



root 




Fig. 5: A graph i9 with p Vo = 2,p Vl = 2,p V2 = 3,p V3 = 2, p V4 = 2 and fc = 12, and some labels. The line 
numbers, distinguishing the lines are not shown. The lines length should be the same but it is drawn of 
arbitrary size. The momentum flowing on the root line is u = v_{vq) = sum of all the nodes momenta, 
including v_ v . 

We attach to each node (or "vertex") v of the tree a vector v_ v e Z d , called a "mode 
label" , and to a line oriented from v to v' we attach a "current" or "momentum" v_ (v) 
and a "'propagator" g{v_(v)): 

K(v)= f E g^o-v^'i^-v)- 2 (4.8) 

if v is the "first node" of the tree, see vq in Fig. 5, then the momentum v_{v) is called 
the total momentum of the tree. 

One defines the value Val ($) of a tree d decorated with the described labels as 

vaiw=«- n f^r^m (4 - 9) 

DSnodes v_ u — y " 

where /,, are the Fourier coefficients of the perturbation f(a) = „ f ve 1 -'— , with 
| y_ | < N, fv = f-v , and v' denotes the node following v and, if v is the first node of 
the tree (so that v' would be the root which (by our conventions) is not a vertex) then 
v_ v i is some unit vector e, see Fig. 3. 

Given the above Feynman rules, which one immediately derives from (4.6), (4.7) the 
component along the vector e, labeling the root, of the k-the order Fourier coefficient 
h ^ of the function h ( ip ) • e , which we write as 

h„ ■ e =eh^ ■ e + e 2 ■§.+..., (4.10) 

is simply the sum of the values Vali? over all trees # which have total momentum y_ , k 
nodes, and no branch carrying -momentum. 

One can check directly that h{ip) so defined is a formal solution of (4.4): the series 
(4.10) with the coefficients defined as above is called the Lindstedt series of the KAM 
problem: it was introduced, at least as a method for computing the low order coefficients 
h^ k \ by Lindstedt and Newcomb in celestial mechanics problems and it was shown to 
be possible to all orders by Poincare (one has to show that the algorithm generating the 
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series does not produce graphs with branches carrying momentum which, by (4.9), 
would yield meaningless expressions for the corresponding tree values), see [Po92], vol. 
3. 

The number of trees of order k that do not differ only by the labeling of the lines (i. e. that 
are topologically the same) is bounded exponentially in k while the total number of trees 
(labels included) is, therefore, of order k\ times an exponential in k. 

Therefore taking into account the in (4.8) we see that the perturbation series might 
have convergence problems only if there exist individual graphs whose value is too large, 
e.g. 0(kV) for some 7 > 0. 

Such graphs do exist; an example: 

—2 

WWW W 2t 

\ \ \ V 

r- I -J 

-HO -HO -Ho -Ho -fc/3 

Fig. 4: A tree of order k with momentum v_ = v l + . . . j£j./ 3 and value of size of order (fc/3)! T if 
u_ = v 1 + . . . vj. is "as small as it can possibly be": or is "almost resonant" i.e. such that w ■ v_ ~ 
C 1 1 I ^~ - The tree has fc/3 + 1 branches carrying momentum v_ = ^^/^ i.e. the k/3+ horizontal 
branches. The last fc/3 — 1 branches have momenta v_ ^, i = 2, . . . , fc/3 so arranged that their sum plus 
u_ 1 , i.e. v_, is "almost resonant". 

Therefore even though the theory is loopless and its perturbation series is well defined 
to all orders, yet it is non trivial because the fc-th order might be "too large" . 

This is a typical situation in "infrared divergences" due to too large propagators: in fact 
the "bad" graphs (like the one in Fig. 5) are such because ( lo ■ v)~ 2 , «- e - the propagator, 
is too large. 

It is remarkable that the same strategy used in the analysis of the Fermi gas theory, 
i.e. the renormalization group approach outlined in general terms in §3, works in this 
case. We decompose the propagator as 

j^Tf = t X ~^f - t W.o ■ *) (4.n, 

h— — oo ' h— — oo 

where if h < we have set x^( x ) = x(2~ hx ) with x(x) = unless x is in the interval 
2 _1 < \x\ < 1 where \( x ) = 1; an d X^ 1 ^ is defined to be identically 1 for \x\ > 1 and 
otherwise so that 1 = J2h=- 00 X^( x ) : see Fig- 2 and (2. 9). 6 

Given a Feynman graph, i.e. a tree z9 (with the decorating labels) one replaces g{ u - v) 
by the last sum in (4.11) and we obtain trees with branches bearing a " scale label' . 



In this problem there is no need to use functions which arc not as in Fig. 4, i.e. with smoothed out 
discontinuities: there is, however, a minor difficulty also in this case. The decomposition above, (4.11), 
can be done exactly as written only if lu () e R d is outside a certain set of zero volume in R d . Although 
already the Diophanntine property (4.3) holds only outside a set of zero volume in R d , as is well known, 
the set of w e R d for which what follows can be done literally as described is slightly smaller (although 
still with a complement of zero volume). However the following discussion can be repeated under the 
only condition that u verifies (4.3) provided one does not insist in taking a sequence of scales that are 
exactly equal to 2 , h = 1, 0, — 1, . . . and one takes a sequence of scales that have bounded but suitably 
variable successive ratios. Here we ignore this problem, see [Ga94] for the cases that work as described 
here (first reference) and for the general case (second reference) discussed under the only "natural" 
condition (4.3). 
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We collect the graph lines into clusters of scale > h with at least one line of scale h and 
define the order n of the cluster to be the number of nodes in it. 

It is not difficult, see [P684] and [Ga94], to see that only graphs which contain one 
incoming and one outgoing line with the same momentum (hence same scale) are the 
source of the problem: for instance in the case of the graph in Fig. 4 all horizontal lines 
have the same momentum of scale h ~ logfc~ r while all the other lines have scale O(l) 
because <n. The subgraphs that contain one incoming and one outgoing line with 
the same momentum of scale h < have been called in [E196] resonances, perhaps not 
very appropriately given the meaning that is usually associated with the word resonance 
but we adopt here the nomenclature of the latter breakthrough work. The subgraphs 
which would be resonances but have maximal scale h = 1 are not considered resonances. 

Therefore it is natural to collect together all graphs which contain chains of 1, 2, 3 . . . 
clusters with equal incoming and outgoing momenta on scale h < 0. For instance the 
graph in Fig. 4 contains a chain of fc/3 clusters, each containing a single line (namely 
the lines with v_ a , — v_Q modes at their extremes), and one cluster with fc/3 — 1 lines (the 
lines entering the node 2fc/3 + 1) and fc/3 + 1 lines external to the resonant clusters, the 
horizontal lines. Call 



v ) = sum of e™ times the value of the clusters with n nodes and single 
-' . 6 (4.12) 

incoming and outgoing lines of equal momentum and scale h 

which can be given meaning because disregarding the incoming line we can regard the 
subcluster with one entering and one exiting line as a tree (or subtree) with root at 
the node where the exiting line ends, so that its value r^ h \u) will be naturally defined 
as the value of $ times (w ■ v_)~ 2 (which takes into account the propagator of the line 
entering the cluster). 

This leads to a rearrangement of the scries for h in which the propagator of a line with 
momentum v_ on scale h is T^ h \v) rather than (uj_ ■ hL)~ 2 X^ h \^.o ' e) and there are 
no more clusters with one incoming and one outgoing line of equal momentum. 

Of course the graphs with k nodes give a contribution to h which is no longer propor- 
tional to e k because they contain quantities T^ h \u) which are (so far formally) power 
series in e. 

The same argument invoked above, [P6] , [Ga94] , gives again that if for some constant 
R and for all e small enough one could suppose that 

\T {h \iL)\<R{^-K)- 2 (4.13) 
then the series for h v would be convergent for e small enough. 

Furthermore if (4.13) holds for scales 0, — 1, . . . , h + 1 then one sees that there exist 
ah,bh, Th such that 

r<*>(,) = + + (1 + r h (^ Uo • ,) (4.14) 

with 1 77, | < Re, \g( h \ui ■ v)\ < N 2 (uj_ ■ v_)~ 2 ', and we see therefore that a/j,6/j play 
the role of "running couplings" : they are non trivial functions of e but they are the only 
quantities to control because if we can show that they are bounded by sR, say, then the 
convergence of the perturbation series would be under control because we are reduced 
essentially to the case in which no graphs with resonances are present. 

Naturally, by the principle of conservation of difficulties, the quatities ah and bh are 
given by power series in e whose convergence is a priori difficult to ascertain. 
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As in the case of the Fermi gas, setting c h — (ah, bh) the very definition of such constants 
in terms of sums of infinitely many Feynman graphs implies that they verify a recursion 
relation 

c h = Ac h+1 +B h {c h+1 ,...,c ), h<0 (4.15) 

where A is a 2 x 2 diagonal matrix with diagional elements 2 2 ,2, and c 1 = because 
by definition ^= only for h < 0. The analogy with the previous Fermi system 
problem seems quite strong. The function B h ({c})h+i....,o * s n °t homogeneous in the 
c h : and this is an important difference with respet to the previous fermionic case. 

Clearly if ah, bh ^ for some h we are "lost" because although, under the boundedness 
condition (4.13), on the running couplings the beta function Bh is well defined by a 
convergent series, and although the whole rearranged perturbation series is convergent 
under the same condition, we shall not be able to prove that the running couplings stay 
small so that (4.13) is selfconsistent. In fact both ah, bh are "relevant couplings" (because 
the elements of the matrix A are > 1) and the two data a\,b\ (which are for h = 1 
because, by definition, there are no resonances of scale 1) will need to be carefully tuned 
so that the renormalization group trajectory ah,bh that (4.15) generates from a\,b\ is 
bounded: there are however no free parameters to adjust in (4-.1)! 

The situation is very similar to the one met in the Fermi liquid theory: in that case one 
solves the problem by showing that the beta function vanishes, at least asymptotically, for 
the marginal couplings A^, Sh- it remains the relevant coupling Vh which can be bounded 
only because we have in that problem the freedom of "adjusting" a free parameter (the 
chemical potential v). 

In the present case we have two relavant parameters, ah,bh, and no free parameter in 
the Hamiltonian: the situation would be hopless unless it just happened that the correct 
initial data for a bounded renormalization group trajectory were precisely the ones that 
we have, namely ai = b\ = 0. This means that one should prove the identity £^(0) = 
for all h < so that c h = Q for all h < 1. 

The vanishing of the beta function at c = was understood in [E196] and it can be 
seen in various ways, see also [Ga94] . It is however always based on symmetry properties 
of the model (the choice of the origin of ip plays the role of a "gauge symmetry" : the 
interpretation of the cancellations used by [E196] as a consequence of a gauge symmetry 
was pointed out and clearly stated in [BGK99], although of course the symmetry is used 
in the analysis in [E196] and in great detail in [Ga94]). 

We are back to a very familiar phenomenon in field theory: a non renormalizable the- 
ory becomes asymptotically free and in fact analytic in the parameter e measuring the 
strength of the perturbation, because of special symmetries which forbid the exponential 
growth of the relevant couplings in absence of free parameters in the Lagrangian of the 
model which could possibly be used to control them. 

Concerning the originality of the results obtained with the techniques exposed in this 
paper the following comments may give an idea of the status of the matter: 

(1) The KAM theory presented above is only a reinterpretation of the original proofs 
[E196], following [Ga94] and [GM96a] see also [BGK98]. However even in classical me- 
chanics the method has generated new results. Other problems of the same type that can 
be naturally interpreted in terms of renormalization group analysis of suitable quantum 
fields, see [GM96a], [GM96b], and in the theory of the separatrix splitting, see [Ga95], 
[GGM99]. 

(2) The results on the theory of the Fermi systems were obtained for the first time by the 
method described above (including the vanishing of the beta function, [BG90]) and have 
led to the understanding of several other problems [BM95], [BGM99], [BM99], [Ma97], 
[Ma98a], [Ma98b], [Ma99]. 
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